function dx=Plant(t,x,flag,para)
dx=zeros(2,1);

ut=para(1);
dt=para(2);


dx(1)=x(2)+dt;
dx(2)=-2*x(1)-x(2)+exp(x(1))+ut;
